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Under the assumption that dark matter is made of new particles, annihilations of those are 
required to reproduce the correct dark matter abundance in the Universe. This process can occur 
in dense regions of our Galaxy such as the Galactic center, dwarf galaxies and other types of sub- 
haloes. High-energy gamma-rays are expected to be produced in dark matter particle collisions and 
could be detected by ground-based Cherenkov telescopes such as HESS, MAGIC and VERITAS. 
The main experimental challenges to get constraints on particle dark matter models are reviewed, 
making explicit the pros and cons that are inherent to this technique, together with the current 
results from running observatories. Main results concerning dark matter searches towards selected 
targets with Cherenkov telescopes are presented. Eventually, a focus is made on a new way to 
perform a search for Galactic sub-halos with such telescopes, based on wide-field surveys, as well as 
future prospects. 

PACS numbers: 95.35.-|-d, 95.85.Pw 

I. INTRODUCTION 

Both from the ground and using satellites, observations of the cosmos in 7-rays have 
known significant improvement in the last twenty years. The Cherenkov technique has now 
reached a mature state, with several running observatories around the globe. The very good 
sensitivities of these instruments allow not only to observe 7-ray sources such as supernovae 
remnants, pulsars or active galactic nuclei, but also to search for exotic signals, in particular 
dark matter induced. In this paper, an attempt is made to describe the minimal knowledge 
that is necessary to have about this experimental technique to understand its pros and cons 
in the quest for non baryonic dark matter. 

In the first section, the technique of Cherenkov telescopes is described, with an emphasis 
on the different methods to reduce backgrounds. The second section is devoted to the 
expectations for possible dark matter signals, the way to choose targets and to display the 
constraints. The third and fourth sections deal with a review of results from all experiments 
using targeted searches and blind searches, as well as some prospects. 

II. GAMMA-RAY ASTRONOMY WITH CHERENKOV TELESCOPES 

A. Principle 

The basic idea of running ground based telescopes to observe cosmic 7-rays is to use the 
atmosphere as a calorimeter. When a high energy particle hits the top of the atmosphere, 
it induces a cascade of secondary particles. At energies of about 1 TeV, that cascade is 
fully contained in the atmosphere, and it produces a flash of Cherenkov photons. Hadronic 
and electromagnetic particles produce different types of cascades. While hadrons induce 
irregular particle showers, electrons, positrons and 7-rays produce a more even shower. The 
Cherenkov flash is contained in a cone of 1° opening angle produced at 10 km height. The 
projection of the cone on the ground is a disk of order 250 m diameter, and the flash 
lasts about 5 ns. From any place inside this disk, the atmospheric shower is observable, 
should one use a sensitive enough instrument. Ground-based 7-ray observatories use this 
principle to measure 7-ray induced Cherenkov light, as sketched on the left panel of Fig. [TJ 
Large dishes are used to collect enough photons, those are focused on very sensitive cameras 
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equipped with photomultipliers. The cameras are able to integrate the signal very quickly 
and to resolve the image of the atmospheric showers, in order to fight against different 
types of backgrounds, as we shall see in the following. To gain in angular resolution, energy 
resolution and background subtraction, several telescopes are used simultaneously to observe 
the event, thus getting a stereoscopic view of the particle cascade. 




FIG. 1: Left: Principle of the Cherenkov telescopes, several telescopes measure the Cherenkov flash 
induced by the initial 7-ray (image from K. Bernohr [T]). Right: Schematic representation of the 
different layers of background that Cherenkov astronomers have to suppress to obtain a image of 
the TeV sources. 



B. Backgrounds, data analysis and caveats of the technique 

As it can be inferred from the description of the technique in the previous subsection, the 
main backgrounds one has to fight against in order to have an image of a TeV source in the 
sky are: 

• The night-sky background : diffuse light from stars, light pollution, the Moon, etc 

• Cherenkov light induced by hadronic cosmic rays 

• Cherenkov light from electron and positrons cosmic rays 

• Diffuse 7-ray emission that could overcome the signal from the source 

Schematically, the image of the source is obtained when all these backgrounds are mas- 
tered, as shown on the sketch of Fig. [T] We shall now describe the basic principles that 
allow to suppress these backgrounds (see for a pedagogical application). 

The night sky background (NSB) is fought against mainly using fast integration elec- 
tronics. The closer the observation to the 5 ns of the actual shower flash, the less NSB 
is integrated. Fig. [2] shows simulations of images of an atmospheric event obtained in the 
camera of a Cherenkov telescope. Here the same simulated 7-ray induced signal is integrated 
over 100 ^s, 1 /iS and 10 ns. For the longest integration times, the NSB signal dominates, 
whereas the electromagnetic shower appears clearly in the camera when the integration is 
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done fast. Depending on the experiments, different techniques are used (ADC flashes, analog 
memories), and the actual integration of the signal is of order 10 ns at best. Another way 
to suppress NSB is to use a triggering strategy that favors signals that are clustered, as the 
NSB is expected to be fairly constant over the field of view of the cameras. 
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FIG. 2: Three simulated image of a 7-ray shower, obtained with different integration times: from 
left to right 100/is, l/is and 10 ns (courtesy K. Bernlohr). 

Other types of background arc mainly suppressed during the offiine analysis. First, one 
has to get rid of the images of cosmic ray induced showers. This is handled by the capability 
of the telescopes to imaging the atmospheric event. Because heavy cosmic rays induce 
irregular showers, their images can be rejected during data analysis. In Fig. [3] different 
images obtained in the focal plane of a telescope are displayed. In this figure, the left panel 
shows an image of a hadronic shower, the central panel is the image of a single muon ring 
(from a hadronic shower that is outside the telescope field of view) , and the right panel is the 
image of a 7-like event. The offiine analysis tools are designed to reject the two first ones, 
and keep the last one, based on the topology of the image. To get an idea of the rejection 
power of the analysis, for a typical bright source, 10 h of observation yield 10'' recorded 
events, out of which 10'' are actual 7-rays. After the topology-based analysis procedure, 
10^ events are selected, dubbed 7-like events. At this point of the analysis, electrons and 
positrons cosmic rays are still present in the sample. 
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FIG. 3: Simulated images of atmospheric events induced from cosmic particles, as observed in the 
focal plane of a Cherenkov telescope. From left to right: hadronic shower, isolated muon, 7-like 
event (courtesy J. Hinton). 
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FIG. 4: Illustration of the ON-OFF background suppression technique (figure from 

At fixed energy, electrons, positrons and 7-rays all induce identical electromagnetic show- 
ers, thus 7-like events contain all these species. A feature of charged TeV scale cosmic 
rays is that their distribution at the top of the atmosphere is isotropic. The subtraction 
technique is then based on the estimate of the isotropic part in a region of the sky where 
no 7-ray signal is expected. This region is called the OFF region, and the corresponding 
event rate is subtracted where the source is expected to be (ON region). Fig. |4] displays two 
possible choices for OFF regions, and a histogram showing the signal from the source and 
the background. The later is composed of electrons, positrons, residuals from misidentified 
hadronic showers. Note that if an actual 7-ray signal is diffuse all over the field of view, 
this technique will likely erase it. The general rule is that a diffuse signals is observable 
with Cherenkov telescopes provided its extension is smaller than the field of view of the 
experiment, which is about a few degrees. Observing larger scales diffuse emission requires 
non standard techniques such as dedicated OFF runs. 

The energy of the initial 7-ray is deduced from the amount of collected Cherenkov pho- 
tons, and comparison to Monte Carlo simulations matching the conditions of data taking. 
Uncertainties lie in the knowledge of the atmospheric conditions, simulation uncertainties 
and intrinsic shower fluctuations. At the end, the energy resolution is about 10%. 

C. Ongoing experiments 

Current generation of Cherenkov telescopes all run in stereoscopic mode. In the northern 
hemisphere, MAGIC uses two large telescopes (17 m diameter) in the Canary islands, and 
VERITAS runs 4 telescopes of 12 m diameter. VERITAS is the successor of the Whipple 
telescope and is located in Arizona. In the southern hemisphere, the Cangaroo collaboration 
runs 4 10 m telescopes in Australia, and HESS consists of 4 telescopes of 12 m diameter in 
Namibia. Fig. [s] shows the locations of the observatories on a map. 
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FIG. 5: Current generation of Cherenkov telescopes observatories around the globe. 

The Cherenkov technique is complementary to satellite telescopes, such as Fermi. Fermi 
works between about 300 MeV and 300 GeV, has a full sky coverage with a large field of 
view, and a relatively small effective area, of order 1 m^. Whereas Cherenkov telescopes 
run only at night when the sky is clear, Fermi can record data all around the clock, that 
compensates for the smaller effective area. The basic features one has to keep in mind in 
order to understand the pros and cons of using Cherenkov telescopes to search for dark 
matter are the following: 

• Large effective areas, of order 10^ m^ 

• Relatively small fields of view, of order a few degrees 

• Angular resolution of about 0.1° 

• Energy Threshold between 50 GeV and 200 GeV 

• 10% energy resolution 

• A duty cycle allowing ~1000 h of data taking a year 

D. The future : CTA 

The next step for 7-ray astronomy will consist of the building of a larger array of Cherenkov 
telescopes -the Cherenkov Telescope Array (CTA) The observatory will take advantage 
of proven experimental techniques on HESS, MAGIC and VERITAS, and use 50 to 100 
telescopes. One site in each hemisphere will allow a full sky coverage. On each site, the 
observatory will be composed of 3 groups of telescopes, from the center of the array to its 
outskirts: large 23 m telescopes, medium 12 m telescopes and small 6 m telescopes, as shown 
on Fig. [6j This is intended to cover a wider energy range from a few 10 GeV to 200 TeV. 
The sensitivity will be improved by a factor of 10 compared to current generation arrays 



FIG. 6: Artist's view of a possible layout for the CTA observatory: 3 classes of telescopes (23 m, 
12 m, 6 m) will allow to cover a wide energy range from 10 GeV to 200 TeV. 

and the angular resolution will be improved by a factor of up to 4. The energy resolution 
and the duty cycle will be similar to what exist now. The field of view will depend on the 
class of the telescope, from ~3° for large ones to up to 8° for small ones. This is a technical 
limitation, as it is difficult to build large telescopes with long focal distance. The layout of 
CTA will allow a flexible operation: deep observations of a region of the sky involving the 
whole array, larger surveys with subarrays pointing in different positions on the sky, regular 
observation with 3/4 of the array while the other quarter monitors variable sources, etc. 

III. PARTICLE DARK MATTER AND ITS INDIRECT DETECTION 
A. Relic density and 2-particles annihilations 

The next half of the paper is devoted to the constraints one can put on particle dark 
matter models thanks to the good sensitivity of Chcrcnkov telescopes. This first section 
recalls the main arguments for particle dark matter and the derivation of possible signals. 

The cosmological standard model stipulates that 84% of the matter in the Universe is 
non-baryonic. This is motivated by different probes on various scales [5]. As an example, 
one can notice that the Cosmic Microwave Background (CMB) homogeneity seems to be in 
conflict with the very existence of galaxies. Indeed, the CMB radiation shows a homogeneous 
Universe at the level of (5 = 5p/ p ~ 10~^. Nowadays obviously 6^1. As (5 grows propor- 
tionally to the scale factor of the metric, if the Universe was as homogeneous as the CMB 
was at recombination, one would have only 5 ^ 10~^ in the present. This simple fact points 
out the need to introduce a new type of matter which could present larger fluctuations. Not 
to affect the CMB, this matter should not interact with light, therefore being non-baryonic. 
In an independent way, models beyond the standard model of particle physics predict the 
existence of new massive stable particles that have the required properties to make up the 
cosmological non-baryonic matter, dubbed Dark Matter (DM) in the following. Interestingly 
enough, these DM candidates were at first not proposed to solve the non-baryonic matter 
issue but appear when trying to address the limitations of the standard model of particle 
physics such as the grand unification or the stabilization of the Higgs mass with respect to 
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radiative corrections. The current cosmological DM density is set by their annihilation rate 
in the early Universe. This provides a natural value for the annihilation cross section of 
{av) ^ 3 X 10^^^ cm^s^^, where {av) is the velocity- weighted annihilation cross-section. 

B. Expected signals in Cherenkov telescopes and derivation of constraints 

In the standard picture, the stars of our Galaxy lie in a thin disk of ^ 20 kpc radius, 
which is dipped into a spherical DM halo of scale a factor of ~ 10 larger. Some of the DM 
particles happen to collide within the halo, producing standard model particles. This exotic 
production of standard model particles is associated with the emission of 7-ray with energies 
of the order of the DM particle mass. The considered weakly interactive massive particle 
(WIMP) being related to electroweak physics, its mass is expected to range between 100 
GeV and a few TeV. The 7-ray emission associated with DM particle annihilations provides 
a chance to detect DM particles through their collisions. 

The expected energy-integrated 7-ray flux from a region of volume V &i & distance D with 
DM density p is 



where is the number of 7-rays per collision, m is the DM particle mass. The number 
is strongly dependent on the energy threshold of the experiment, as it is the integral 
of all photons from its threshold to the DM mass. As the expected DM induced photons 
energy spectra are quickly falling with energy, it is then always better to have a low energy 
threshold, in order to be sensitive to more photons. The defined volume V can contain a 
whole DM structure {e.g. a in dwarf galaxy) or part of it (like for the Milky- Way halo). 

Up to now, no convincing DM signal has been observed, so the way to present the results 
is to draw a contour in the av-m plane. The flux sensitivity of an experiment is determined 
by the characteristics of the corresponding observation: effective area, exposure time, data 
quality, zenith angle of the target etc. The experimental work consists precisely in deter- 
mining the flux sensitivity. Then, for a fixed particle DM mass m, a limit on the cross 
section av is computed for the expected signal not to be larger than the sensitivity. This 
can be performed provided the DM density within the observed region is properly modeled. 
This is the main source of uncertainty in this type of analysis. The constraints are usually 
compared to model predictions that are computed with codes such as DarkSUSY [B] or 
micrOMEGAs [7] . Depending on the particle physics model, the dependence of iV^ on the 
mass of the DM particle is different. Some collaborations (Whipple, VERITAS, sometimes 
HESS) choose to use a generic parameterization for N^im) which is not exact but somehow 
gives the average value over a possible range of models. It is also possible to extract the 
value of iV^ from the model and compute the constraint point by point. That is the way 
MAGIC presented its results and in that case the constraint is not a continuous line. An- 
other approach is to show the constraints obtained with the two extreme final sates which 
are r+r^ and bh. 

As for the choice of the targets, the haloes of galaxies and clusters of galaxies are assembled 
through the merging of a huge number of smaller structures. Most mergers are incomplete 
and large CDM halos, e.g. the one around the Milky Way, harbor an enormous population 
of subhaloes, which are a record of its assembly history. Some of these subhaloes, for the 
most massive of them, contain baryonic matter and stars, making up the satellites of the 
Milky- Way. A careful study of the kinematics of those stars allow obtaining constraints 
on the DM mass content of the satellite, which can help to reduce the uncertainty on the 
constraint. However, the presence of baryons in the center of haloes can modify the density 
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profile of tfie DM. Then, only subhaloes that do not contain baryons are expected to have 
their primordial shape untouched. The interaction of the satellite or the subhalo with 
the gravitational potential of the Milky- Way can also cause the profile to be significantly 
modified by tidal effects. In this review we shall see results from all types of sources: dwarf 
galaxies with strong tidal disruption or none of it, globular clusters for which baryons could 
have played a role, the Galactic center where there is an actual conventional 7-ray source, 
and DM clumps, which are basically untouched subhaloes. In all cases, the ideal signal one 
would expect from a dense DM region has the following features: 

• Slightly extended, with an extension of order 1° and a halo-type morphology 

• Steadiness: it should have strictly no time variability 

• No high energy counterpart 

• No nearby conventional source 

• Cutoff at electroweak-like scale in the energy spectrum 

• A few alter-egos with the exact same properties in the sky 

C. Searches towards known targets 

The most used strategy is to select targets where the WIMP annihilation process can 
occur efficiently, and observe them with a sufficient amount of time. Targets are selected 
thanks to the observation of stars and study of their kinematics. This allows to infer the 
DM mass content of the target, although with usually large uncertainties. The knowledge of 
their positions in the sky is very precious for Cherenkov telescopes, because of their relatively 
small fields of view and background subtraction techniques. The presence of baryons can 
have drawbacks too, they can alter the DM content of the target or induce conventional 
7-ray emissions. 

The most dense region of our Galaxy is its center, located at 8.5 kpc. HESS dedicated 
almost 50h of observation to the Galactic center 8j. Among all the targets discussed here, 
it is the only one for which a 7-ray emission has actually been found. It has been shown 
that most of the emission cannot be due to DM annihilations. The reason is that the energy 
spectrum of the source does not fit a DM-like spectrum. One should note also that unlike 
the other targets mentioned here, several astrophysical objects that could be at the origin of 
the emission are indeed present in the same field of view. To infer constraints in the {av}-m. 
plane, the authors of |8] search for the maximal cross section for which the energy spectrum 
does not exclude the corresponding DM contribution. Results are displayed in Fig. |9j in 
which the band represents the uncertainty from the particle physics models on Nj. 

Other observations have been conducted by HESS, Veritas, MAGIC and Whipple. All 
of them get their constraints from the fact that no signal has been found. Fig. [7] presents 
a summary some of the observed targets, with their types, distances, the corresponding 
experiment and exposure used to set the constraint. Only objects for which constraints 
are drawn in Fig. |8] and Fig. [9] appear in the figure, other objects have been used for 
searching for particle DM, such as Bootes [12], M33 and M32 fT3| and the Perseus galaxy 
cluster jT^, but those are not listed here. The results from different analyses are shown 
on Fig. |8] and Fig. |9] They are sort into two groups, depending on wether the DM content 
modeling suffer from important systematic uncertainties or not. The first group, on Fig. [8] 
concerns targets that are not strongly affected by tidal effects. Fig. [9] displays constrains from 
somehow more optimistic analyses. It includes observations of the Canis Major overdensity. 
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FIG. 7: Galactic coordinates representation of the targets for which the constraints are displayed in 
Fig.jsjand Fig. [9] and their main characteristics: distance, type and experiments that observed them, 
with corresponding exposure. The displayed targets are the Galactic center (HESS [F), Sagittarius 
dwarf (HESS 9 ), Sculptor (HESS [13]), Carina (HESS ,10 ), Draco (MAGIC [TT], VERITAS [H], 
Whipple 113]), Willman I (VERITAS T?, MAGIC [14]), Ursa Minor (VERITAS TS', Whipple [13]), 
Canis Major (HESS ^15^) and M15 (Whipple [S]) 



Sagittarius dwarf galaxy (that experienced strong tidal stripping), the results from the 
Galactic center, as well as the results from M15, for which a very optimistic assumption has 
been made (adiabatic contraction around a hypothetical central black hole) . As explained in 
the previous section, the constraints have different shapes because of the parameterization 
of and the different energy dependence of the effective area of the experiments. For some 
of the results, an uncertainty band represents the variation of the constraints when the DM 
content of the target is varied within the region that is compatible with star kinematics 
(called halo uncertainty in the legends). For the sake of comparison, the result from Fermi 
on Sculptor dwarf galaxy is shown on Fig. |8] One can note that the constraints are very 
complementary. The advantage of Fermi here is the much lower energy threshold (below 
1 GeV) and the very long exposure. In the case of "safe" dwarves, the limits on (crw) lie 
in the 10"^"'^ — 10"^"^ cm^s~^ region. For the more uncertain targets the limits lie between 
10~^^ cm^s"^ and 10~^^ cm^s^^ 

The next generation of Cherenkov telescopes, CTA, will certainly conduct the same type 
of observation. The effective area will be a factor of 10 larger, and constraints will benefit 
from a lower threshold. In il7^ though, the authors show that only for optimistic particle 
physics scenarios CTA will be sensitive to natural WlMPs. 
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HESS Carina (incl. halo uncertainty) 
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FIG. 8: Compilation of constraints in the av~m plane, from HESS, Veritas, MAGIC, Whipple, in 
the case of targets poorly affected by tidal stripping (Carina, Draco, Sculptor, Ursa Minor, Wilmann 

I)- 
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HESS Galactic center (incl. model uncertainty) 



Whipple IVI15+adiabatic compression 



HESS Canis Major (incl. halo uncertainty) 
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FIG. 9: Compilation of constraints in the av-m plane, from HESS, Veritas, Whipple, in the case 
of targets for which more uncertainties are associated (Canis Major overdensity. Galactic center, 
M15+adiabatic contraction, Sagittarius dwarf galaxy) 



D. Blind searches for DM clumps 

Although Cherenkov telescopes cannot perform wide- field surveys in a single shot, the 
maturity of this technique allows now to conduct scans of significant fractions of the sky, 
and make blind searches for new sources. This has been done with the HESS experiment, 
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FIG. 10; Constraints obtained from the search for DM clumps with HESS wide field galactic plane 
survey (left panel) and prospects for CTA (right panel) in the case of a galactic plane survey and 
a survey of one fourth of the sky. 



for scanning the Galactic plane [TH]. This type of scan allows searching for DM clumps. 
Indeed the Milky Way halo is the result of the merging of a large number of smaller haloes, 
some of which are still present in the Milky Way, as predicted in high resolution N-body 
simulations [12]. In ^0], data from the HESS survey has been used to build a sensitivity 
map. The authors of [ST] have shown that a large number of subhaloes are predicted to be 
in the scanned region (displayed on Fig. [t] as a dashed rectangle) . For given values of av 
and m, annihilation within those objects could make a statistically significant number of 
them shine enough in 7-rays to be observable. From the observation perspective, none of the 
unidentified sources discovered in the HESS scan present the required characteristics to be 
DM clumps. Then, from the convolution of the sensitivity map with the prediction from DM 
clustering in numerical simulations, it is possible to get the constraints that are shown on 



the left panel of Fig. 10 As shown on this figure, the results are again very complementary 
to Fermi wide field searches [22] although Fermi has a full sky coverage. The drawback 
from a much smaller scanned regions is compensated in the case of Cherenkov telescopes 
by a better flux sensitivity. The obtained constraints are of the same order of magnitude 
as for targeted searches. Note however that the systematic errors on the DM distribution 
-although still present- are very different in the two cases. 

Prospects for blind searches for DM subhaloes take advantage of the fact that large surveys 
will be conducted by CTA independently of the search for DM. At least 2 large scans can 
be foreseen: a Galactic plane survey and a wider survey of one fourth of the sky. In [5T| . 
constraints on DM models are computed assuming no DM subhalo candidate is found. In 
the case of a HESS-like survey, constraints are computed with the same exposure. That 
would allow for an improvement of the HESS constraints by a factor of 10 (see right panel 
of Fig. 10 1. To go further, the 1/4 sky survey will have to be used. In that case, in [21] 
the scan regions has been optimized: centered on the Galactic center, with the exclusion of 
latitudes between ±1.5°, (blue region in Fig.[7|. With the same assumptions, and within 6 
years of data taking, this allows to reach the thermally produced WIMPs region, as shown 



on right panel of Fig. 10 
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IV. CONCLUSIONS 



Cherenkov telescopes now allow for deep searches for particle dark matter both from iden- 
tified targets and performing blind searches. Those two types of methods have very different 
types of systematic errors and reach the same level of constraints. CTA will certainly give 
the opportunity to dig deeper into the models, in particular in the light of the future results 
from current observations by the Fermi satellite. 
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